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Two-dimensional electron systems at oxide interfaces are often influenced by a Rashba type spin- 
orbit coupling, which is tunable by a transverse electric field. Ferromagnetism near the interface 
can simultaneously induce strong local magnetic fields. This combination of spin-orbit coupling and 
magnetism leads to asymmetric two-sheeted Fermi surfaces, on which either intra- or inter-band 
pairing is favored. The superconducting order parameters are derived within a microscopic pairing 
model realizing both, the Bardeen-Cooper-Schrieffer superconductor with inter-band pairing, and a 
mixed-parity state with finite-momentum intra-band pairing. We present a phase diagram for the 
superconducting groundstates and analyze the density of states, the spectra, and the momentum 
distribution functions of the different phases. The results are discussed in the context of super- 
conductivity and ferromagnetism at LaA103-SrTi03 interfaces and superconductors with broken 
inversion symmetry. 
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I. INTRODUCTION 

Reducing the symmetry of a superconducting system 
changes the properties of the possible superconducting 
(SC) states. In particular, if the inversion symmetry 
is broken, the combination of Rashba spin-orbit cou- 
pling (SOC) and magnetic fields leads to new classes 
of SC states. Superconductivity with SOC was first 
discussed in inversion-symmetry breaking heavy fermion 
systems 1 ' 2 , and this topic recently expanded in the con- 
text of topological superconductors 3 and in normal-SC 
heterostructures 4 ' 5 . With the discovery of supercon- 
ductivity at the interface between LaA103 (LAO) and 
SrTiOs (STO) 6 ' 7 , a system was identified which mani- 
festly exhibits a combination of SOC, magnetism, and 
superconductivity. At this interface, a dilute, almost 
ideally two-dimensional electron liquid forms 8 10 . The 
polarity of the interface induces a perpendicular electric 
field E and gives rise to a Rashba SOC 7(E x p) • S, 
where p and S are momentum and spin operators, re- 
spectively. Quantitative estimates for the SOC param- 
eter 7 are under debate 11 15 . Recent experiments also 
revealed inhomogeneous in-plane magnetism at the in- 
terface 16 18 which coexists with superconductivity at low 
temperatures and is possibly generated through oxygen 
vacancies 19 ' 20 . A physically similar situation may also be 
realized in thin films of heavy- fermion superconductors in 
an external magnetic field. It is this special concurrence 
of superconductivity with magnetism and SOC which al- 
lows for the formation of unanticipated multi-component 
SC states. 

The magnetic field B couples through a Zeeman term 
/ieB • S to the spin in the same way as E x p through 
the Rashba SOC. Both couplings result in a band split- 
ting and a two-sheeted Fermi surface. Various studies 
focused on the emergence of topological edge states in 



two-dimensional systems with B || E 21 ' 22 . However, in- 
teresting new physics emerges in a boundary free system 
with an in-plane field B _L E. An in-plane field does not 
generate vortices in the superconducting state, but cou- 
ples to the electrons through a Zeeman term. While for 
B = intra-band pairing with zero center-of-mass mo- 
mentum (COMM) is expected, a finite COMM is required 
if the in-plane component of the magnetic-field exceeds 
a certain limit 23 ' 24 . If the magnetic field dominates over 
SOC, a crossover to inter-band pairing must occur, of ei- 
ther BCS type with zero COMM for low fields, or a finite 
COMM state above a critical field, as proposed by Fulde 
and Ferrell 25 and by Larkin and Ovchinnikov 26 . The 
regime with both SOC and magnetic field has previously 
been treated within a Ginzburg-Landau analysis for the 
"helical phase" of the intra-band pairing regime 27 , which 
was later extended to other non inversion-symmetric sys- 
tems 28 ' 29 . Surface superconductivity in magnetic fields 
was discussed also on the basis of the linearized Gor'kov 
equations 23 ' 30 . 

In this paper we characterize the possible SC ground- 
states in the presence of both SOC and in-plane magnetic 
field within a microscopic model. We describe the tran- 
sition from the intra-band pairing regime with dominant 
SOC to the inter-band regime with a dominant Zeeman 
coupling. These two regimes are typically separated by 
a first-order phase boundary and exhibit distinctly dif- 
ferent signatures in the density of states (DOS). On the 
basis of existing experimental data for the LAO-STO in- 
terface we discuss the implications for the SC state which 
may be obtained from tunneling spectra. 

II. THE NORMAL STATE 

For the microscopic description of a two-dimensional 
electron system we use a single-band tight-binding model 
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FIG. 1: (Color online) Fermi surfaces with SOC and magnetic 
field for the tight-binding dispersion ek and band filling n — 
0.24. (a) a = 0.4 1, H x = 0. (b) a = 0.4 1, H x = 0.3 1. (c) 
a = 0, H x = 0.4t. 

on a square lattice with TV sites. We include a nearest- 
neighbor hopping t and a next-nearest-neighbor hopping 
t' with a kinetic energy of the form Tio = s e^c^c^g 
with ek = — 2t(cos k x + cos k y ) + 4£' cos k x cos k y — fi; n is 
the chemical potential. 

The spin part Hs of the Hamiltonian consists of 
the Rashba SOC and the Zeeman coupling to an in- 
plane magnetic field. With the spin operator = 
J2 S s , c\ s cr ss >Ci S > /2 on the lattice site z, Hs becomes in 
momentum space 

n s = ^2(A Ks +H X + isH y ) 4, a Qk,- fl (1) 

with ylk,s = tt(smk y + iss'mk x ), a = ft 2 |E|7/2, H = 
/ieB, and s = d=l. The effective electron mass is denoted 
by m* . For free electrons moving in a perpendicular elec- 
tric field, the parameter 7 = e/(2m* 2 c 2 ) derives directly 
from the Dirac Hamiltonian. However, the effective value 
of 7 can be strongly enhanced in multi-band systems with 
atomic SOC 14 ' 15 . 

The term Hs generates superpositions of electronic 
states with spin t and | and equal momentum k, rep- 
resented by the fermionic operators a\a = + </>kCk| 
and <2k2 = — </>kCk| with the phase factor 

= A^ + H x + iH y 
^ \A^ + H x +iH y [ y) 

A momentum dependent in-plane spin orientation re- 
sults. The indices 1 and 2 denote the two energy bands 
£ki and £k2 resulting from the non-interaction Hamilto- 
nian, which is diagonal in the operators <2ki and a^' 

Ko + = ^2 (^kl«kl a kl + ^k2« k 2 a k2) (3) 
k 

with f k i >2 = e k ± \ A\ef + H x + iH y \. 

The two corresponding Fermi surfaces are shown in 
Fig. 1 for a charge density n = 0.24. Three qualita- 
tively different cases are identified: (a) For finite SOC 
(ft ^ 0) but no magnetic field (H = 0), two concentric 
Fermi surfaces emerge, on which the spins are oriented 
anti-parallel (£ki) or parallel (^2) to E x Vek, i.e. the 



spin winds once around its quantization axis upon circu- 
lating the Fermi surface, (c) In a magnetic field (here H 
in x-direction) with vanishing SOC (ft — 0), the Fermi 
surfaces are also concentric, however with spins oriented 
anti-parallel (£ki) or parallel (^2) to H. In the follow- 
ing we shall mostly analyze the intermediate regime with 
both « / and H x ^ shown in (b). The Fermi sur- 
faces are displaced from the Brillouin-zone center per- 
pendicular to H, e.g. to the positions Ki = (0,Ki) and 
K2 = (0, K2) for H y = 0. The offsets are maximal on 
the critical line a = H X1 and in general \K\\ > |i^2|- 
For small magnetic field s H x <C ft , the offsets are given 
by Ki = K 2 ~ H x /^a 2 + 4tj/xj in first order of H x 
(c.f. Refs. [23,24]). The spin- winding number switches 
from zero for a < H x to ±1 for ft > H x . A first-order 
transition between two superconducting states of differ- 
ent symmetry class is therefore expected at this crossover. 



III. THE PAIRING INTERACTION 

In case (a), intra-band pairing of anti-parallel spins 
with zero COMM is expected on both Fermi surfaces. 
Although spin-triplet inter-band pairing with a finite 
COMM is allowed in the presence of SOC, it is unlikely 
to be energetically favorable. In (c) only singlet pairing 
is possible for a local attraction. This requires inter- 
band pairing and allows for the formation of a Larkin- 
Ovchinnikov (LO) state with finite COMMs q and — q, 
where q connects the two Fermi surfaces, if |H| is larger 
than the SC order parameter 41 . Which pairing type 
is realized in the intermediate regime (b) with a 7^ 
and H x ^ is difficult to anticipate. Nevertheless, if 
H x <C ft we expect an intra-band pairing state where fi- 
nite COMMs qi and q2 are realized which optimize pair- 
ing on band 1 and band 2, respectively. These COMMs 
should be chosen such that the momenta of the paired 
eigenstates lie on the Fermi surface of the corresponding 
band, i.e. qi « 2Ki and q2 ~ 2K 2 . (The exact values of 
qi and q2 can only be found from a self-consistent numer- 
ical calculation and may deviate slightly from the expec- 
tations formulated above.) Such a state was considered 
in Ref. [29], where it was assumed that only one COMM 
exists which optimizes pairing on the larger Fermi sur- 
face. 

In the following we analyze the pairing originating 
from an on-site pairing interaction and characterize the 
crossover from intra- to inter-band pairing for increasing 
field H x . The interaction term in momentum space has 
the form 

k,k' q s 

Expressing Hi in terms of the band operators <2k,i and 
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fflk,2 leads to 
V 



Hi 



2N 2 



[Kven(k, k', q) + V odd (k, k', q)] 



k,k',q 

a,P,a',P' 



a L a -k+q/J a -k'+q/3'ak'a' 



(5) 



where 7/3/3/ = 1 for f3 = f3' and 7/3/3/ = —1 for j3 ^ /?'. 

In the operators <2ki and ak2, the interaction is non- 
local and acquires the momentum dependence of (j) k 
(with p-wave symmetry for H = 0). The interac- 
tion is decomposed into a component V r eve n(k, k', q) = 
Fg* ven (k, q)^ even (k / , q) which is even under the permu- 
tation k o — k + q and a component lodd(k, k', q) = 
Vg* dd (k, q)g dd(k', q) which is odd under the same per- 
mutation, where 



# even (k,q) 
#odd(k,q) 



-k+q) /2, 
-k+q) /2. 



(6) 
(7) 



The mean-field pairing order parameters (OPs) 
A a /3(k, q) generated by the interaction (5) must 
be antisymmetric under an exchange of the paired 
quasiparticles, i.e. 

A a/3 (k,q) = -A /3a (-k + q,q). (8) 

From this condition follows that the two intra-band OPs 
are of odd parity: 



An(k,q) 
A 22 (k, q) 



k' 

V 



N 



V r odd(k, k', q)(a_k'+ q ia k 'i>, (9) 



^ Kdd(k, k', q)(a_ k ' + q2«k'2), (10) 



whereas inter-band pairing allows for odd and even parity 
OPs: 

A? 2 dd (k,q) = A° dd (k,q) 

= ^^^ odd ( k ' k/ ' q ^ a - k/ + ( i iak/2 )' ( n ) 
k' 

A- e "(k,q) = -A|r n (k,q) 



x 21 

-E 

k' 



l/even(k,k / ,q)(a_k'+qiak'2). (12) 



In the mean-field decoupling of Eq. (5), all OPs coupling 
to the intra-band term ^k'c^-k'+qa: mus ^ be of odd parity 
with the OP Af2 6n (k, q) canceling out. The inter-band 

term ^'c^-k'+q/? (a ^ (3) allows for both, odd and even 
parity OPs. However, all odd parity OPs cancel out in 
these terms. The mean-field pairing interaction therefore 
takes the form 



^ F = ^E[ A ° dd ( k ^) 

l (k,q) (ai.al 



k,q 



+A e 



k+ql 



k+q2 



a k2 a -k+q2 , 



a k2 a -k+ql 



(13) 



with 

A odd (k,q) = [A n (k,q) - A 22 (k,q)]/2, (14) 
A eve "(k,q) = A- en (k,q). (15) 

These OPs may be expressed through the k-independent 
quantities 

A eve„/odd (q) = 5even/odd(k) q) A CTe "/° dd (k, q), (16) 

since ^even/odd (k, q) cancels the k-dependence in 
A eve„/odd (k)q) K 5 * ven/odd (k,q). In the limit H x 
one finds A„(k,q) = -A 22 (k,q) and A even (k,q) = 0, 
thus a pure odd-parity OP is realized. 

The OPs {A 11 (k,q),A 22 (k,q),A? dd (k,q)} form a 
triplet in the band space spanned by aki and <2k2, and 
Af2 6n (k) is a band-space singlet. They can be trans- 
formed back to the standard spin singlet -0(k, q) and 
triplet d(k, q) = (^(k, q), d y (k, q), d z (k, q)) OPs 1,31 us- 



ing the inverse transformation 
one finds that 



^(k,q) = g, 



even 
even 



(k,q)A eve "(k,q) 
(q) + A odd (q) 



Ck s . In this way 

- 5odd (k,q)A odd (k,q) 
(17) 



and 

4(k,q) 



5odd (k,q)A even (k )q )+ 5e 



1 (k,q)A odd (k,q), 
(18) 

whereas the equal-spin OPs <ic(k, q) and c^(k, q) are ab- 
sent for the local pairing interaction. Singlet pairing is 
thus realized if either g even (k, q) or g dd(k, q) is zero, i.e., 
if either a or H x is zero. In this case, the parity of the OP 
is either purely even or purely odd. If both, a and H x 
are finite, the parity is mixed and therefore a spin-triplet 
component is admixed. 

A phase-coherent superconducting state is realized if 
one (or a few) COMMs are macroscopically occupied. 
We choose a test set of COMMs according to the expec- 
tations on pairing we formulated above, i.e., two COMMs 
qi and q2 should be present to optimize intra-band pair- 
ing in band 1 and band 2, and two additional COMMs 
±q which are required to form an LO state. Further, 
the COMM q = is likely to represent the groundstate 
for small values of a and H x and is therefore included as 
well. Possibly existing higher harmonics of these vectors 
are neglected here. 

In the presence of Q different COMMs q^ in the test 
set, the Hamiltonian % m = H + 7-L$ + %Y * s expressed 
by a 2(Q + 1) x 2(Q + 1) matrix as 



7T = J2 A l 



\Ak, qQ 



Ak, qi 

-^-k+qi 





^qg 



\ 



A k 



-£-k+q Q / 



(19) 
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intra-band 


a 


H x 


qi 


q2 


An(qi) 


A 22 (qi) 


An(q 2 ) 


A 22 (q 2 ) 


al 


0.4 


0.1 


(0,0) 


(0,0) 


-0.234 


0.297 


-0.234 


0.297 


a2 


0.4 


0.223 


(0,0) 


(0,0) 


-0.132 


0.178 


-0.132 


0.178 


bl 


0.4 


0.25 


(0,-21) 


(0, 18) 


-0.060 


0.044 


-0.048 


0.109 


b2 


0.4 


0.3 


(0, -28) 


(0,27) 


0.002 


0.001 


-0.012 


-0.035 


inter-band 


a 


H x 


qi 


q2 


A^ en (qi) 


A? 2 dd (qi) 


A^ en (q 2 ) 


A?f (q 2 ) 


c 


0.1 


0.2 


(0,0) 


(0,0) 


0.225 




0.225 




d 


0.05 


0.44 


(-52,0) 


(52,0) 


0.080 




0.080 





TABLE I: Summary of solutions from the four distinct SC phases for V = 1.7 1 (except 2-band b: V = 2.4 1), t' — 0, and 
n — 0.24, calculated on a 640x640 lattice. For intra-band pairing, An(qi, 2 ) and A 22 (qi )2 ) correspond to pair-momenta qi, 2 , 
whereas for inter-band pairing, as well as for the LO solution, only the OPs Ai 2 en (qi j2 ) exist, a, H and all OPs are given in 
units of t and qi )2 in units of 27r/640a. 



0.4 
0.2 



intra-band 
al 

(q = 0) 



a2. 



normal 
metal 



inter-band (q = 0) 




HJt 0.2 



0.4 



these results and compare them to the above expectations 
in terms of an a-H x phase diagram (Fig. 2). 

In the presence of at least two different COMMs, 
the system is inhomogeneous with finite off-diagonal el- 
ements of the density matrix 34 . In the presence of the 
two COMMs qi and q2, it exhibits a charge-stripe or- 
der with wavevector qi — q 2 . Within the mean- field de- 
coupling of the interaction (5), additional order param- 
eters, representing charge order, could be defined. How- 
ever, their amplitude remains small for the values of the 
COMMs relevant here, and therefore we subsequently ne- 
glect these order parameters. 



IV. RESULTS 



FIG. 2: (Color online) Qualitative phase diagram for the 
model defined by the Hamiltonian (13) for V — 1.7 1 and 
t' — 0.2 1, and n — 0.24 1. The blue regions in the upper- 
left half indicate intra-band pairing with dominating SOC, 
whereas the pink regions in the lower-right half represent 
inter-band pairing with either q x = or q x ^ (LO). Darker 
color in direction parallel to the a — H x line indicates an 
increasing pair momentum q y . The stars mark the positions 
of the five examples listed in Table I. The critical magnetic 
field if L o = A even (H x = 0) marks the breakdown field of the 
classical BCS superconductor. 



With A£ = (a k l,ak2,a t _ k+qil ,aI k+qi2 ,...,aI k+qQl , 

«L k+qQ2 )- Here, E^ and A kq are the 2x2 matrices 



and 



*k,q 



£k 



A odd (k,q) 
-A even (k,q) 




A even (k,q) 
-A odd (k,q) 



(20) 



(21) 



has to be diagonalized numerically to calculate the 
OPs A even / odd (q^) self-consistently. Below we analyze 



We choose the parameters n = 0.24 and V = 1.7 1, 
for which stable solutions for all the different SC phases 
exist. Figure 2 shows the phase diagram in the a- 
H x plane. As expected, intra-band pairing dominates 
throughout the half-plane a > H x: although supercon- 
ductivity is steadily suppressed by an increasing mag- 
netic field. If H x > 0, the normal-state Fermi sur- 
faces shift as shown in Fig. 1 and two eigenstates on 
opposite points of one Fermi surface differ in energy 
by AE(k) « |f a (k F ) - £ a (k F - K a )\. Nevertheless, 
if AE(k) < A odd (q = 0) for all k, pairing with zero 
COMM is favored (phase a). In lowest order in H X1 one 
finds that the maximal value of A£?(k) (for k x = 0) is 
2aH x /s/a 2 + 4t|/x|. If this value exceeds A odd (0), only 
solutions with finite COMM pairing exist (phase b). This 
phase is similar to the "helical phase" of Ref. [29], al- 
though in our calculations two COMMs qi ~ 2K X and 
q 2 ~ 2K2 coexist in the groundstate of this phase, op- 
timizing pairing in band 1 and band 2, respectively (so- 
lution bl in Tab. I). Upon approaching the boundary of 
phase b, pairing on the smaller Fermi surface (of band 1) 
is suppressed (solution b2) and a one-band SC state re- 
mains. 42 The finite pair density in band 2 induces an en- 
ergy gap also in band 1 . This "interior gap" phenomenon 
was predicted also in two-band systems with strongly dif- 
fering effective electron masses 32 ; in our model however, 



-0.6n 0.6/r ~-0.6n 0.6/r 




-0.6tt 0.6tt -0.6tt 0.6/r 



FIG. 3: (Color online) Gap structure of the superconducting solutions of Table I. For each k x , the eigenenergies Ek x (k y ) are 
plotted as a function of k y (except for the LO solution (d), where Ek y (k x ) is plotted). The eigenenergies of band 1 and band 2 
are plotted in pink and in blue, respectively. The opacity of each point encodes the weight with which the corresponding state 
contributes to the density of states, given through the coherence factors. 



the induced energy gap is momentum dependent (see dis- 
cussion of Fig. 3 below). 

Inter-band pairing is possible below a = H x . The 
phase transition at a = H x is first order and typically 
ends at a critical point, above which intra- and inter-band 
pairing are separated by a non-superconducting metallic 
phase. The position of this critical point changes with 
system parameters such as V, t' and the electron density, 
and occurs at a = H x = 0.082 1 in Fig. 2. At the critical 
magnetic field H LO (a) (with H LO (0) = A even (q = 0)), 



a first-order transition into the LO phase takes place, 
which is stable within a certain field range above i^LO 
and well separated from intra-band pairing. The extent 
of the SC phases on the H x axis scales with the interac- 
tion strength V. The qualitative structure of the phase 
diagram remains however unchanged upon changing V. 

In the following we discuss the characteristic proper- 
ties of the two types of intra- and the inter-band pairings 
on the basis of four typical parameter sets, marked with 
white stars in Fig. 2 and summarized in Table I. As 
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nil 



P (k) 




FIG. 4: n(k) (left panel) and p(k) (right panel) for the so- 
lutions al) - d) of Tab. I. For the arrows in bl) and d) see 
main text. 



already mentioned above, superconducting states with 
at least two coexisting COMMs are generally inhomo- 
geneous, i.e., states with COMMs qi and q2 form a 
"pair-density wave" with wave vector (qi — q2)/2 and 
a concomitant charge-density wave with a wave vector 
qi — q2 34 ' 36 - The properties of the different phases 
are essentially determined by the structure of the mo- 
mentum dependent energy gap resulting from the OPs 
A even (° dd )(k,q). This is analyzed in Fig. 3, showing the 
eigenenergies Ek x (k y ) for all k x as a function of k y (ex- 
cept Fig. 3 (d) showing E^ y (k x )). Figure 3 (al) shows the 
zero COMM intra-band pairing state with an almost uni- 
form energy gap A odd (k, 0) in both bands. Because the 
normal-state Fermi surfaces are shifted in /^-direction, 
the "gaps" are tilted as a function of k y . As a conse- 
quence, the coherence peak in the density of states (DOS) 
splits into a double-peak structure. 

Upon increasing a or H x , the displacement of the 
Fermi spheres and therefore the tilting of the energy gaps 
increases, reaching a point where the upper and lower 
bands touch [Fig. 3 (a2)]. Beyond this limit, zero COMM 
superconductivity breaks down and is replaced in a first- 
order transition by a finite COMM state [Fig. 3 (bl)]. In 
this state, pairing with finite COMM q 1? 2 ~ 2Ki 5 2 com- 
pensates the shift of band 1 or of band 2, respectively, 
and generates an energy gap covering the full Fermi sur- 
face of band 1 or of band 2. On the other hand, these 
pairings also induce an interior gap in the other band, 
which results in an overall energy cost, because the inte- 
rior gap shifts more eigenstates to higher energies than to 
lower energies. Therefore, upon reaching the outer phase 
boundary of the SC phase b, pairing on the smaller Fermi 
surface (band 1) is lost, leaving only the gap in band 2 
and the induced interior gap in band 1 [Fig. 3 (b2)]. 

The gap structure of the zero COMM inter-band pair- 
ing state is shown in Fig. 3 (c). Clearly visible is a gap 
shifted to higher energy in band 1 and a gap shifted to 
lower energy in band 2. At the critical magnetic field i^ c , 
the gap edges in band 1 and band 2 touch, thus for mag- 
netic fields above H C1 zero COMM pairing is no longer 
possible. The LO state above H c [Fig. 3 (d)] has a gap 
structure similar to the finite COMM intra-band pairing 
state. However, because of the required inter-band pair- 
ing, there is no COMM generating an energy gap which 
covers the whole Fermi surface of one band. Instead, the 
COMMs ±q of the LO state are equal to the momen- 
tum difference of the Fermi surfaces of band 1 and band 
2 (c.f. Fig. 4 (d), where paired states are indicated by 
the two-colored arrows). Each band therefore contains 
two tilted energy gaps for +q and — q, which overlap 
the Fermi energy only partially. The direction of the 
COMMs of the LO state is not a priori clear. While the 
states with q in x- and in ^/-direction are equivalent for 
a = 0, in the calculations with a = 0.05 1 (see Table I), 
COMMs parallel to H are favored. 



A complementary characterization of the distinct SC 
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phases is given by the momentum distribution 

n(k) = (4 f c kt ) + (c^c H ) (22) 
(left panels of Fig. 4) and the pair density 

p 2 (k) = ^|(c_ k+qi c kt )| 2 (23) 
q 

(right panels of Fig. 4). Fermi surfaces are represented by 
lines of discontinuities in n(k), which becomes continuous 
wherever the Fermi energy falls into an energy gap. The 
pair density p(k) is typically largest where the derivative 
of n(k) is largest (these are the dark red and dark blue 
regions in n(k)), but it is zero where n(k) changes dis- 
continuously. Pairing occurs only where an energy gap 
crosses the Fermi energy, but no pairing occurs at interior 
gaps, since these do not influence n(k). Interior gaps are 
induced by pairing on the Fermi surface of the other band 
and do not contribute to the condensation energy. In the 
intra-band pairing states (al, a2, bl, and b2), pairing oc- 
curs between two blue quasi-particle states and between 
two red quasi-particle states, whereas in the inter-band 
pairing states [c and d] , pairing occurs between a red and 
a blue quasi-particle state (indicated by the two-colored 
arrows in Fig. 4 d). 

The low-energy physics is characterized essentially by 
the OP A even /° dd (q). Most eminent is the fact that n(k) 
and p(k) are almost x-y symmetric in both zero COMM 
pairing phases (al and a2), although x-y symmetry is 
clearly broken in the normal state. On the other hand, 
the finite COMM SC states (bl, b2, and d) break x- 
y symmetry. This change of symmetry upon increasing 
the magnetic field might be observable experimentally, as 
we will argue below. 

We have thus identified four distinct SC phases, the 
different symmetry properties of which are most evident 
in their corresponding pair densities p(k) (right panels 
of Fig. 4). Characteristic for the finite COMM phases 
is also a fascinating fine structure in both n(k) and 
p(k), originating from partially reconstructed Fermi sur- 
faces (marked by black arrows in Fig. 4 bl) and Fermi 
pockets (blue arrows in Fig. 4 bl), which are typically 
present in an inhomogeneous SC system with two differ- 
ent COMMs 34 36 . These structures are most likely not 
observable experimentally. A characteristic quantity in 
which they may leave fingerprints is the DOS (Fig. 6). 
For the chosen band filling (n = 0.24), the normal- state 
DOS is featureless near the Fermi energy. Zero COMM 
pairing typically generates a full energy gap, although a 
double-peak structure in the coherence peaks appears in 
a finite magnetic field. This double peak originates from 
the tilting of the dispersion. At the critical magnetic field 
i^LCb the inner coherence peaks touch and the gap closes. 
Characteristic for finite COMM pairing is a finite DOS at 
the Fermi energy and additional structure away from the 
Fermi energy which originates from the interior gaps. In 
the LO state, parts of the Fermi surface around k x = 
do not lie within the energy gap and therefore survive 
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FIG. 5: (Color online) Density of states as a function of the 
energy relative to the Fermi level for dominating SOC (upper 
panel) and dominating magnetic field (lower panel). Upper 
panel: intra-band solutions al) (dark blue), bl) (light blue), 
and b2) (dashed purple) of Table I. Lower panel: inter-band 
solution c) (light blue) and LO solution d) (dark blue) of 
Table I. 



into the SC state. The origin of the in-gap states in the 
finite COMM intra-band pairing phase is more complex; 
they originate from the partially reconstructed Fermi sur- 
faces. These in-gap states are visible in Fig. 3 (bl) close 
to k y = and are responsible for the lines of disconti- 
nuities crossing the dark red and dark blue "circles" in 
Fig. 4 (bl) (marked by the black arrows) 43 . 



V. DISCUSSION 

A prominent example, in which the combination of su- 
perconductivity with SOC and magnetism is relevant, is 
the superconducting LAO-STO interface. We therefore 
estimate which SC phases may occur in this system. The 
parameters characterizing the LAO-STO samples vary 
considerably in different experiments. Here we use the 
values given by Caviglia et al. for the density of mobile 
electrons at the interface and the Rashba SOC 7 ' 11 . The 
electron density for zero gate voltage is 4.5 x 10 13 cm -2 , 
which corresponds to n ~ 0.07 mobile electrons per unit 
cell. For a quadratic dispersion relation one obtains 
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FIG. 6: (Color online) Sketch of the phase diagram for the ex- 
perimentally obtained parameters of the LAO-STO interface 
system (see main text). Here, B x — H x /{ib and a e x P = oq, 
where a is the lattice constant of STO. For the smallest re- 
ported value of ttexp (10 -12 eVm), we estimate the transition 
into the finite COMM intra-band pairing phase to occur at 
~ 5 T, while a — H x corresponds to ~ 40 T. A reentrance into 
the inter-band pairing phase upon increasing H x is therefore 
not possible. 

E F « 35meV and the energy scale t « 80meV. The 
measured values of the Rashba SOC strength a exp = aa 
(where a is the lattice constant of LAO) vary with the 
gate voltage from 1 to 5xl0 -12 eVm 11 . In our model, 
these values correspond to a « 0.03 - 0.15 1. 

The intrinsic magnetism at the interface is strongly in- 
homogeneous, thus no definite position on the H x axis 
can be given. The average in-plane magnetic moment 
varies between 0.01/iB 18 and 0.3/xb per unit cell 17 . This 
corresponds to ~ 15 - 500 G, or to H x w 0.1 - 3 x 10" 5 t 
in our model. For all of these values, the system is 
therefore deep in the zero COMM intra-band phase. For 
a^> H x , the phase bound ary to the f inite COMM pairing 
phase occurs at H x w A^a 2 + 4t\/n\/2a w 0.05-0.25^, 
where A is the SC energy gap. From measurements of the 
superfluid density 38 , an estimate A « 40ueV can be de- 
duced. The in-plane magnetic field required for the finite 
COMM phase (phase b) is therefore around 1-5 T, which 
is far above the measured internal magnetic field, but can 
easily be reached in the laboratory. The LAO-STO in- 
terface is therefore a likely candidate for the observation 
of a finite COMM superconducting state. 

By applying an even larger external magnetic field one 
possibly might wonder whether this system can be moved 
into the inter-band phase. The magnetic field H x = a 
corresponds to ~ 40 T for a = 0.03 1. According to our 
discussion above, this is the minimal field above which 
this inter-band pairing is possible. However, this field 
is far larger than the critical field H^q = A even (H x = 
0) « 40 ueV. For the LAO-STO interface, the energy gap 
is therefore ~ 10 2 - 10 3 times smaller than the spin- 



orbit band splitting reported in Ref. [11] (and the order 
parameters obtained in Table I) and a reentrance into the 
inter-band superconducting state is not possible for this 
system. The phase diagram obtained from our model for 
the experimentally determined parameters of the LAO- 
STO interface system is thus quantitatively different from 
the one shown in Fig. 2. In particular, superconductivity 
extends only little in the direction of the magnetic field 
and the normal state is reached for magnetic fields H x <C 
a. This situation is schematically shown in Fig. 6. 

Theoretical estimates for a are difficult. A simple 
electrostatic model for the interface, using Dime's result 
a = \E\eti 2 1 '4m* 2 : c 2 yields a value for the Rashba SOC 
strength that is many orders of magnitude smaller than 
the values reported from the experiments. In this sce- 
nario, inter-band pairing is expected even in moderate 
external magnetic fields. However, for electron densities 
above ~ 1.5-1.9 x 10 13 cm -2 , a multi-band behavior has 
been inferred from Hall effect measurements 14 , for which 
the effective value of a can be strongly enhanced 14 ' 15 . 
By tuning the electron density at the interface through 
a gate voltage 9 , both the one-band and the multi-band 
regime may be accessible 14,37 . Therefore, superconduc- 
tivity with weak and strong SOC may both be realized at 
the LAO-STO interface, depending on the gate-voltage 
tuned electron density. 

Because of the uncertainties in the parameter values 
and the sensitivity of our model to the strength of the 
magnetic field, it is difficult to predict details in the 
structure of the DOS. For the reported internal magnetic 
fields, a double-peak structure is unlikely to be resolvable 
experimentally, and was indeed not seen in Ref. [38]. 

In some heavy- fermion superconductors like CePtsSi 
or CeIrSi3 inversion- symmetry is intrinsically broken 2 . 
The physics discussed above may emerge in thin films of 
these materials in an external magnetic field, but is al- 
tered by vortices in the presence of a perpendicular field 
component 29 . In the case of CePtsSi, the SOC band 
splitting is of the order of 100 me V 39 and exceeds the 
size of the SC energy gap by several orders of magnitude. 
This system would realize therefore the intra-band pair- 
ing phase and again, no reentrance into the inter-band 
pairing phase is possible. A correct description of the 
heavy-fermion superconductors is however more compli- 
cated, since their pairing interaction is non-local. This 
allows also for equal-spin triplet pairing 40 and will there- 
fore lead to a far richer phase diagram which has yet to 
be explored. 
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